SHOCK WAVES AND COMPACTONS FOR 
FIFTH-ORDER NONLINEAR DISPERSION EQUATIONS. II 



VICTOR A. GALAKTIONOV 

Abstract. The following first problem is posed: to justify that the standing shock wave 

S-(x) = — signx = {—1 for x < 0, 1 for x > 0}, 

is a correct "entropy" solution of the Cauchy problem for the fifth-order degenerate 
nonlinear dispersion equations, NDEs (same as for the classic Euler one u t + uu x = 0) 

u t = —{uu x ) xxxx and u t = —{uu xxxx ) x in R x R + . 

These two quasilinear degenerate PDEs are chosen as typical representatives, so other 
(2m + l)th-order NDEs of non-divergent form admit such shocks waves. As a related 
second problem, the opposite initial shock S+ (x) = — S-(x) = sign x is shown to be a non- 
entropy solution creating a rarefaction wave, which becomes C°° for any t > 0. Formation 
of shocks leads to nonuniqueness of any "entropy solutions" . Similar phenomena are 
studied for a fifth-order in time NDE Uuttt — {uu x ) xxxx in normal form. Other NDEs, 

e.g., u t = -{\u\u x ) xxxx + \u\u x in Ixl + , 

are shown to admit smooth compactons, as oscillatory travelling wave solutions with 
compact support. The well known nonnegative compactons, which appeared in various 
applications (first examples by Day, 1998, and Rosenau-Levy, 1999), are nonexistent in 
general and are not robust relative small perturbations of parameters of the PDE. 

This is more extended and detailed version of the arXiv preprint |19j . Particularly, 
essential novelties are available in § 5.4, where a family of similarity extensions after 
blow-up was detected, which were mentioned but not found in [19j § 5]. 



1. Introduction: nonlinear dispersion PDEs, and main directions of study 

1.1. Five main problems and layout: shocks, rarefaction waves, and com- 
pactons for fifth-order NDEs. Let us introduce our basic models, which are five fifth- 
order nonlinear dispersion equations (NDEs). These are ordered by numbers of derivatives 
inside and outside the quadratic differential operators involved on the right-hand sides: 



(1.1) 


u t = 




(NDE-(5,0)) 


(1.2) 


u t = 


(uu xxxx ^j x 


(NDE-(4,1)), 


(1.3) 


u t = 


(^uu xxx j xx 


(NDE-(3,2)), 


(1.4) 


u t = 


(uu xx ) xxx 


(NDE-(2,3)), 


(1.5) 


u t = 


(uu x ^j xxxx 


(NDE-(1,4)). 
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The only fully divergent operator is in the last NDE-(1,4) that, being written as 

(1.6) u t = -{uu x ) xxxx = -\ {u 2 ) xxxxx (NDE-(1, 4) = NDE-(0, 5)) , 

becomes also the NDE-(0,5), or simply the NDE-5. This completes the list of such 
quasilinear degenerate PDEs under consideration. 

The main feature of these degenerate odd-order PDEs is that they admit shock and 
rarefaction waves, similarly to the first-order conservation laws such as Euler's equation 

(1.7) u t + uu x = in E x R + , u(x,0) = u (x) in E. 

Before explaining the physical significance of the NDEs and their role in general PDE 
theory, we pose four main problems for the above NDEs (the same as for ( 11. 71) ): 

(I) Problem "Blow-up to S- " (Section EJ): to show that the shock of the shape — signx can 
be obtained by blow-up limit from a smooth self-similar solution U-(x,t) of (II. ip - (II. 51) in 
E x (0,T), i.e., the following holds: 

(1.8) u4x,t)^S4x) = -signx={j 1 f °J r X x <°^ as t -> T~ m Lj oc (R). 

(II) The Riemann Problem S + (RP+) (Section [3]): to show that the initial shock 

— 1 for x < 0, 
1 for x > 0, 

for NDEs (I T. jp ~ (j 1.5\) generates a "rarefaction wave", which is C°° -smooth for t > 0. 

(III) The Riemann Problem S_ (RP-) (SectionHJ): introducing a u 5-entropy test" (smooth- 
ing of discontinuous solutions at shocks via a u 5 - deformation" ) , to show that 

(1.10) S-{ x ) i s an "5 -entropy" shock wave, and S+(x) is not. 

(IV) Problem: nonuniqueness / entropy (Section [5]): to show that a single point "gradient 
catastrophe" for the NDE (I i.5j) leads to the principal nonuniqueness of a shock wave exten- 
sion after singularity. This also suggests nonexistence of any proper entropy mechanism 
for choosing any "right" solution after single point blow-up. 

In Section [61 we discuss these problems in application to other NDEs including the 
following rather unusual one: 

Ututt — (uu x ) xxxx , 

which indeed can be reduced to a first-order system that, nevertheless, is not hyperbolic, so 
that modern advanced theory of ID hyperbolic systems (see e.g., Bressan [1J or Dafermos 
[TU] ) does not apply. The main convenient mathematical feature of (11.111) is that it is 
in the normal form, so it obeys the Cauchy-Kovalevskaya theorem that guarantees local 
existence of a unique analytic solution and makes easier application of our 5-entropy 
(smoothing) test. Regardless this, (II. lip is shown to create in finite time shocks of the 
type S-(x) in (II .8p and rarefaction waves for other discontinuous data ~ S+(x) in (II. 9p . 

Finally, we consider the last: 
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(V) Problem "Oscillatory Smooth Compactons" (Section [7]): to show that the perturbed 
version of the NDE ( \1.5h , as a typical example, 

(1.12) u t = -{\u\u x ) xxxx + \u\u x in R x R + , 

admits compactly supported travelling wave (TW) solutions of changing sign near finite 
interfaces. Equation fll.l2f) is written for solutions with infinitely many sign changes, by 
replacing u 2 by the monotone function \u\u. 

Nonnegative compact structures have been known since beginning of the 1990s as com- 
pactons (Rosenau-Hyman, 1993, [Sj). We show that more standard in literature non- 
negative compactons of fifth-order NDEs such as (I1.12p are nonexistent in general, and, 
moreover, these are not robust (not "structurally stable"), i.e., do not exhibit continu- 
ous dependence upon the parameters of PDEs (say, arbitrarily small perturbations of 
nonlinearities). 

1.2. A link to classic entropy shocks for conservation laws. Indeed, the above 
problems (I) (III) are classic for entropy theory of ID conservation laws from the 1950s. 
It is well recognized that shock waves first appeared in gas dynamics that led to mathemat- 
ical theory of entropy solutions of the first-order conservation laws and Euler's equation 
fll.7p as a key representative. The entropy theory for PDEs such as (\1.7\i , with arbitrary 
measurable initial data uo, was created by Oleinik [361 EI] and Kruzhkov [31] (analogous 
scalar equations in ~R N ) in the 1950-60s; see details on the history, main results, and 
modern developments in the well-known monographs [H [TUJ, [50] . Note that first analysis 
of the formation of shocks for (11.71) was performed by Riemann in 1858 |42j; see further 
details and the history in [3]. It is worth mentioning that the implicitly given solution 
u = u(x,t) of the Cauchy problem (11.71) . via the characteristic formula 

u = u (x — ut), 

containing the key wave "overturning" effect, was obtained earlier by Poisson in 1808 |39j ; 
see [30]. 

According to entropy theory for conservation laws such as (11.71) . it is well-known that 
(11.101) holds. This means that 

(1.13) u_(x,t) = S-(x) = — signx 

is the unique entropy solution of the PDE (II. 7p with the same initial data S-(x). On 
the contrary, taking S^-type initial data (11.91) in the Cauchy problem (II. 7p yields the 
continuous rarefaction wave with a simple similarity piece-wise linear structure, 

{—1 for x < —t, 
~ for \x\ < t, 
1 for x > t. 

Our first goal is to justify the same conclusions for the fifth-order NDEs, where, of 
course, the rarefaction wave in the RP+ is supposed to be different from that in (11.141) . 

We now return to main applications of the NDEs. 
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1.3. NDEs from theory of integrable PDEs and water waves. Talking about odd- 
order PDEs under consideration, these naturally appear in classic theory of integrable 
PDEs from shallow water applications, beginning with the KdV equation, 

(1.15) u t + uu x = u XXX) 
the fifth- order KdV equation, 

Uf -\- u xxxxx ~\~ 30 u u x -\- 20 u x u xx -\- 10 uu xxx = 0, 

and others. These are semilinear dispersion equations, which being endowed with smooth 
semigroups (groups), generate smooth flows, so discontinuous weak solutions are unlikely, 
though strong oscillatory behaviour of solutions is typical; see references in [24"| Ch. 4]. 

The situation is changed for the quasilinear case. In particular, consider the quasilinear 
Harry Dym equation 

(1.16) u t = u 3 u xxx , 

which is one of the most exotic integrable soliton equations; see [24"l § 4.7] for survey 
and references therein. Here, (11.161) indeed belongs to the NDE family, though it seems 
proper semigroups of its discontinuous solutions (if any) have never been examined. On 
the other hand, moving blow-up singularities and other types of complex singularities of 
the modified Harry Dym equation, 

_ 3 13 
Ut U U xxx U x 2 U , 

have been described in [6\ by delicate asymptotic expansion techniques. 

In addition, integrable equation theory produced various hierarchies of quasilinear 
higher-order NDEs, such as the fifth-order Kawamoto equation [30], as a typical example 

(1.17) Ut u u xxxxx -\- 5 u u x u xxxx -\- 10 m u xx u xxx . 

We can enlarge this list talking about possible quasilinear extensions of the integrable 
Lax's seventh-order KdV equation 

u t + [35u 4 + 70(u 2 u xx + u(u x ) 2 ) + 7(2uu xxxx + 3(u xx ) 2 + Au x u xxx ) + u xxxxxx ] x = 0, 

and the seventh-order Sawada-Kotara equation 

u t + [63u 4 + 63(2u 2 u xx + u(u x ) 2 ) + 1\[uu xxxx + [u xx ) 2 + u x u xxx ) + u xxxxxx ] x = 0; 

see references in [2U p. 234]. 

The modern mathematical theory of odd-order quasilinear PDEs is partially originated 
and continues to be strongly connected with the class of integrable equations. Special 
advantages of integrability by using the inverse scattering transform method, Lax pairs, 
Liouville transformations, and other explicit algebraic manipulations have made it possible 
to create a rather complete theory for some of these difficult quasilinear PDEs. Nowadays, 
well-developed theory and most of rigorous results on existence, uniqueness, and various 
singularity and non-differentiability properties are associated with NDE-type integrable 
models such as Fuchssteiner-Fokas-Camassa-Holm (FFCH) equation 

(1.18) (I - D 2 x )u t = -3uu x + 2u x u xx + uu xxx = -(I - D 2 x ){uu x ) - [u 2 + \ (u x ) 2 ] . 
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Equation (11.181) is an asymptotic model describing the wave dynamics at the free surface 
of fluids under gravity. It is derived from Euler equations for inviscid fluids under the 
long wave asymptotics of shallow water behaviour (where the function u is the height of 
the water above a flat bottom). Applying to ( 11.18P the integral operator (J — D 2 ) -1 with 
the L 2 -kernel uo(s) = le - ' 5 ' > 0, reduces it, for a class of solutions, to the conservation 
law fll.7p with a compact first-order perturbation, 

(1.19) u t + uu x = -[u* {u 2 + \{u x ) 2 )] x . 

Almost all mathematical results (including entropy inequalities and Oleinik's condition 
(E)) have been obtained by using this integral representation of the FFCH equation; see 
the long list of references given in [24, p. 232]. 

There is another integrable PDE from the family with third-order quadratic operators, 

(1.20) u t - u xxt = auu x + (3u x u xx + uu xxx (a, (3 e R), 

where a = — 3 and (3 = 2 yields the FFCH equation fll . 18j) . This is the Degasperis-Procesi 
(DP) equation for another choice a = —4 and (3 = 3: 

(1.21) u t - u xxt = -4uu x + 3u x u xx + uu XXX} or u t + uu x = - [u * (§ u 2 )] x . 

On existence, uniqueness (of entropy solutions in L 1 H BV), parabolic e-regularization, 
Oleinik's entropy estimate, and generalized PDEs, see [5]. 

Note that, since the non-local term in the DP equation ( 11.211) does not contain u x , the 
differential properties of its solutions are distinct from those for the FFCH one (11. 19ft . 
Namely, the solutions are less regular, and (11.211) admits shock waves, e.g., of the form 

^shock(M) = ~\ signs e~\ x \ 

with rather standard (induced by (11.71) ) but more involved entropy theory; see [321 [14] . 

Besides (11. 18ft and (11 .2111 . the family ( 11.201) does not contain other integrable entries. 
A list of more applied papers related to various NDEs is also available in [2H Ch. 4]. 

1.4. NDEs from compacton theory. Other important applications of odd-order PDEs 
are associated with compacton phenomena for more general non-integrable models. For 
instance, the Rosenau-Hyman (RH) equation 

(1.22) ut = (u 2 ) xxx + (u 2 ) x 

has special important applications as a widely used model of the effects of nonlinear 
dispersion in the pattern formation in liquid drops [47]. It is the K(2, 2) equation from 
the general K(m, n) family of the following NDEs: 

(1.23) ut = (u n ) xxx + (u m ) x (u>0), 

that describe phenomena of compact pattern formation, [431 144] . Such PDEs also appear 
in curve motion and shortening flows [46]. Similar to well-known parabolic models of the 
porous medium type, the K(m, n) equation (11.231) with n > 1 is degenerate at u = 0, 
and therefore may exhibit finite speed of propagation and admit solutions with finite 
interfaces. The crucial advantage of the RH equation (11.22)1 is that it possesses explicit 

5 



moving compactly supported soliton-type solutions, called compactons [13113], which are 
travelling wave (TW) solutions to be discussed for the PDEs under consideration. 

Various families of quasilinear third-order KdV-type equations can be found in [3], 
where further references concerning such PDEs and their exact solutions are given. Higher- 
order generalized KdV equations are of increasing interest; see e.g., the quintic KdV 
equation in [28], and also [55], where the seventh-order PDEs are studied. 

More general B(m, k) equations, 

u t + a(u m ) x = n{u k ) xxx , 

which coincide with the K(m, k) after scaling, also admit simple semi-compacton solutions 
[38] . The same is true for the Kq(m,uj) nonlinear dispersion equation (another nonlinear 
extension of the KdV) [43J 

u < + K) I + K- w (n\)J I = 0. 
Setting m = 2 and uj — ~ yields a typical quadratic PDE 

(1.24) B(u) = u t + {u 2 ) x + uu xxx + 2u x u xx = 0. 

It is curious that (11.241) admits an extended compacton-like dynamics on a standard 
trigonometric-exponential subspaces, on which 

(1.25) u(x, t) = Co(t) + Ci(t) cos Ax + C^t) sin Ax G W$ = Spanjl, cos Ax, sin Ax}, 



where A = 




. This subspace is invariant under the quadratic operator B in the usual 



sense that B(W 3 ) C W 3 . Therefore substituting (fl~25|) into the PDE 03} yields for the 
expansion coefficients on W3 {Co, C\, C2} a 3D nonlinear dynamical system; see further 
such examples of exact solutions of NDEs on invariant subspaces in [2"4"| Ch. 4] . 

Combining the K(m,n) and B(m, k) equations gives the dispersive-dissipativity entity 
DD(k,m,n) [15] 

Ut + a(u m ) x + (u n ) xxx = u(u k ) xx , 

which can also admit solutions on invariant subspaces for some values of parameters. 
For the fifth-order NDEs, such as 

(1.26) u t = a{u 2 ) xxxxx + (3{u 2 ) xxx + j{u 2 ) x inlx R + , 

compacton solutions were first constructed in [UJ, where the more general K(m,n,p) 
family of PDEs, 

u t + Pi(u m ), + Hu n ) xxx + foDl{u p ) = 0, 

with m, n, p > 1, was introduced. Some of these equations will be treated later on. Equa- 
tion (11.261) is also associated with the family Q(l,m,n) of more general quintic evolution 
PDEs with nonlinear dispersion, 

(1.27) u t + a(u m+1 ) x + u[u{u n ) xx ] x + 5[u{u l ) xxxx ] x = 0, 

possessing multi-hump, compact solitary solutions [4~9] . 
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Concerning higher-order in time quasilinear PDEs, let us mention a generalization of 
the combined dissipative double- dispersive (CDDD) equation (see, e.g., [JT]) 

(1-28) u tt = au xxxx + f3u xxtt + j(u 2 ) xxxxt + 5{u 2 ) xxt + e(u 2 ) t , 

and also the nonlinear modified dispersive Klein-Gordon equation (mKG(l, n, k)), 

(1.29) u tt + a(u n ) xx + b{u k ) xxxx = 0, n, k > 1 (u > 0); 

see some exact TW solutions in [29J. For b > 0, (11.291) is of hyperbolic (or Boussinesq) type 
in the class of nonnegative solutions. We also mention related 2D dispersive Boussinesq 
equations denoted by B(m,n, k,p) |54j, 

(u m ) tt + a(u n ) xx + (3{u k ) xxxx + l(u p ) yyyy = in K 2 xi 

See [2H Ch. 4-6] for more references and examples of exact solutions on invariant subspaces 
of NDEs of various types and orders. 

1.5. On canonical third-order NDEs. Until recently, quite a little was known about 
proper mathematics concerning discontinuous solutions, rarefaction waves, and "entropy- 
like" approaches, even for the simplest third-order NDEs such as (11.221) or (see [23| [T8]) 

(1.30) u t = {uu x ) xx . 

However, the smoothing results for sufficiently regular solutions of linear and nonlinear 
third-order PDEs are well know from the 1980-90s. For instance, infinite smoothing results 
were proved in [7j (see also |27j ) for the general linear equation 

(1.31) u t + a(x,t)u xxx = (a(x,t) > c> 0), 
and in [5] for the corresponding fully nonlinear PDE 

(1-32) u t + f{u xxx ,u xx ,u x ,u,x,t) = ( f Uxxx > c> ) ; 

see also [2] for semilinear equations. Namely, for a class of such equations, it is shown 
that, for data with minimal regularity and sufficient (say, exponential) decay at infinity, 
there exists a unique solution u(x,t) £ for small t > 0. Similar smoothing local in 
time results for unique solutions are available for equations in IR 2 , 

(1.33) u t + f(D 3 u, D 2 u, Du, u, x, y, t) = 0; 

see [33] and further references therein. 

These smoothing results have been used in [18] for developing a kind of a 5-entropy 
test for discontinuous solutions by using techniques of smooth deformations. We will 
follow these ideas applied now to shock and compacton solutions of higher-order NDEs 
and others. 

2. (I) Problem "Blow-up": EXISTENCE OF SHOCK SIMILARITY SOLUTIONS 

We now show that Problem (I) on blowing up to the shock S-(x) can be solved in 
a unified manner by constructing self-similar solutions. As often happens in nonlinear 
evolution PDEs, the refined structure of such bounded and generic shocks is described in 
a scaling-invariant manner. 
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2.1. Finite time blow-up formation of the shock wave S~(x). One can see that all 
five NDEs fll.ip - fll.5l) admit the following similarity substitution: 

(2.1) u.{x,t) = g{z), z = x/(-t)* (t<0), 

where, by translation, the blow-up time in reduces to T = 0. Substituting (12. II) into the 
NDEs yields for g the following ODEs in M, respectively: 



(2.2) 




(2.3) 


{gg^)' = -\g'z 


(2.4) 


(g 9 "')" = -lg'z, 


(2.5) 


{gg'T = -\9'z, 


(2.6) 


{ggT = -\9'z 



with the following conditions at infinity for the shocks S-\ 
(2.7) 0(=foo) = ±1. 

In view of the symmetry of the ODEs, 



(2i 



9 >-> -g, 

Z I — > —z. 



it suffices to get odd solutions for z < posing anti-symmetry conditions at the origin, 
(2.9) £(0)=/(0)=s (4) (0) = 0. 

2.2. Shock similarity profiles exist and are unique: numerical results. Before 
performing a rigorous approach to Problem (I), it is convenient and inspiring to check 
whether the shock similarity profiles g(z) announced in (12.1 p actually exist and are unique 
for each of the ODEs fl2.2l) - (l2.6p . This is done by numerical methods that supply us 
with positive and convincing conclusions. Moreover, these numerics clarify some crucial 
properties of profiles, which will determine the actual strategy of rigorous study. 

A typical structure of this shock similarity profile g(z) satisfying (12. 2p . (12. 9p is shown 
in Figure [TJ As a key feature, we observe a highly oscillatory behaviour of g(z) about ±1 
as z — > =Foo, that can essentially affect the metric of the announced convergence (11. 8p . 
Therefore, we will need to describe this oscillatory behaviour in detail. In Figure [2], we 
show the same profile g(z) for smaller z. It is crucial that, in all numerical experiments, 
we obtained the same profile that indicates that it is the unique solution of (12. 2p . (12. 9p . 

Figure [3](a)-(d) show the shock similarity profiles for the rest of NDEs fll.2p - fll.5p . They 
differ from each other rather slightly. 

Remark: on regularization in numerical methods. For the fifth-order NDEs, this 
and further numerical constructions are performed by MatLab by using the bvp4c solver. 
Typically, we take the relative and absolute tolerances 



(2.10) 



Tols = 10~ 4 . 
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Shock similarity profile for NDE-(5,0): u=-u u 

3 r I xxxxx 



gg (5) =-g'z/5 




-200 -150 -100 -50 




Z 



50 100 150 200 



FIGURE 1. The shock similarity profile g(z) as the unique solution of the prob- 
lem ([22D, CUD; z e [-200,200]. 



Shock similarity profile for NDE-(5,0): u=- u u 

t xxxxx 




FIGURE 2. The shock similarity profile g(z) as the unique solution of the prob- 
lem flZ2]), (JZS]); z £ [-20,20]. 



Instead of the degenerate ODE (12.21) (or others), we solve the regularized equation 



(2.11) 
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(5) 



. (\g'z). 



with the regularization parameter v — 10 



where the choice of small v is coherent with the tolerances in (12.101) . Sometimes, we will 
need to use the enhanced parameters Tols = v = 10~ 7 or even ~ 10~ 9 . 
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Shock similarity profile for NDE-(4,1 ): u t =-(u Shock similarity profile for NDE-(3,2): u t =- (u u xxx ). 




(c) equation (|2.5p (d) equation (|2.6p 

FIGURE 3. Shock similarity profiles as solutions of (|2.3p - (|2.6p . (|2.9p respectively. 
For comparison, dotted lines denote the profile from Figures Q] and [2j 



2.3. Justification of oscillatory behaviour about equilibria ±1 and other asymp- 
totics. Thus, the shock profiles g(z) are oscillatory about +1 as z — > — oo. In order to 
describe these oscillations in detail, we linearize all the ODEs (I2.2l) - (l2.6p about the regular 
equilibrium g(z) = 1 by setting g = 1 + g to get the linear ODE 

(2.12) B* s g = -gW-±g'z = 0. 

Note that this equation reminds us of that for the rescaled kernel F(z) of the funda- 
mental solution of the corresponding linear dispersion equation, 

(2.13) Ut = —u xxxxx in Ix R + . 

The fundamental solution of the corresponding linear operator + D\ in (I2.13P has the 
standard similarity form 

(2.14) b(x,t) =t~*F(y), with y = x/t 1 ^ 5 , 
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where F(y) is a unique solution of the ODE problem 

(2.15) B 5 F = -F^ + \ (Fy)' = in R, / F = 1, or = \ Fy on integration. 

However, the operator B 5 in (12.151) is not identical to that in (12.121) . Moreover, this Bg is 
adjoint to B 5 in some indefinite metric and both the operators possess countable families 
of eigenfunctions, which particularly are generalized Hermite polynomials for Bg. We will 
not use this Hermitian spectral theory later on, so refer to [T7J § 9] and [201 § 8-2] for 
further results and applications. 

Let us return to the linearized ODE (12.121) . Looking for possible asymptotics as z — > 
— oo yields the following exponential ones with the characteristic equation: 

(2.16) g{z) ~ ^ => a 4 = $. 

Finally, choosing the purely imaginary root of the algebraic equation in (12.161) with Re a = 
gives a refined WKBJ-type asymptotics of solutions of (12.21) : 

(2.17) g(z) = 1 + \z\~s [Asm(ao\z\i) + B cos(a |z|3)] + ... as z — > — oo, 

where A and B are some real constants satisfying A 2 + B 2 ^ 0. 

The asymptotic behaviour (I2.17P implies two important conclusions: 

Proposition 2.1. The shock wave profiles g(z) solving (jff.jfl - ()£. 61) . flff.Tp satisfy: (i) 

(2.18) g(z) — 1 ^ L 1 (M_), and 

(ii) the total variation of g(z) (and hence ofu_(x,t) for any t < 0) is infinite. 

Proof. Setting |z|f = v in the integrals below yields by ( 12. 171) : 

OO 5 oo 

(i) I \9(z) -l\dz ^ J i^flldz ^ J l -^0-dv = oo; and 

2.19 

v ' +00 00 5 00 

(ii) |#(-)|Tot.Var. = / \g'(z)\dz ~ / dz ~ / ^0- dv = oo. □ 

— oo 

This is in striking contrast with the case of conservation laws ( II. 7p . where finite total 
variation approaches and Helly's second theorem (compact embedding of sets of bounded 
functions of bounded total variations into L°°) used to be key; see Oleinik's pioneering 
approach [36]. In view of the presented properties of the similarity shock profile g(z), the 
convergence in (jl.8p takes place for any i6R, uniformly in R\ (— //, /i), /i > small, and 
in L^ oc (R) for p G [1, oo), that, for convenience, we fix in the following: 

Proposition 2.2. For the shock similarity profile g(z) the convergence (1 1 . S\) with T = 0: 

(i) does not hold in L 1 (R), and 

(ii) does hold in Lj oc (M) ; and moreover, for any fixed finite I > 0, 

(2.20) ||«_(. > t)-S , _(.)|Ui(-/ I i ) = O((-t)B)-^0 as t - 0~. 

n 



Proof of (12.201) is the same as in (12.191) with a finite interval of integration: for I = 1, 



(-t)-vs Mr 1/4 

1 5 5 1 T 1 

II ' Hz, 1 (—1,1) ~ ( — t) 7, J | cos Z4 1 dz ~ (— t) 5 j v~ io I cosf I dt> ~ (— t) s . □ 

Finally, note that each ^(z) has a regular asymptotic expansion near the origin. For 
instance, for the first ODE (12.21) . there exist solutions such that 

(2.21) g(z) = Cz + Dz 3 -^- Q z 5 + ^z 7 + ..., 

where C < and D e K are some constants. The local uniqueness of such asymptotics 
is traced out by using Banach's Contraction Principle applied to the equivalent integral 
equation in the metric of C(— /z, /i), with fi > small. Moreover, it can be shown that 
(12.211) is the expansion of an analytic function. Other ODEs admit similar local represen- 
tations of solutions. 

We now need the following scaling invariance of the ODEs (I2.2p - (I2.6I) : if gi(z) is a 
solution, then 

(2.22) 9a( z ) — a5 9i{a) * s a solution for any a / 0. 

2.4. Existence of a shock similarity profile. Using the asymptotics derived above, 
we now in a position to prove the following: 



Proposition 2.3. The problem < \2.7\) , (\2.9h for ODEs <\2.2X\ - <\2.6} admits a solution g{ 
which is an odd analytic function. 



z 



Uniqueness for such higher-order ODEs is a more difficult problem, which is not studied 
here, though it has been seen numerically. Moreover, there are some analogous results. We 
refer to the paper [25] (to be used later on), where uniqueness of a fourth-order semilinear 
ODE was established by an improved shooting argument. 

Notice another difficult aspect of the problem. Figures [TH3] above, which were obtained 
by careful numerics, clearly convince that the positivity holds: 

(2.23) g(z) > for z < 0, 

which is also difficult to prove rigorously; see further comments below. Actually, (12.231) 
is not that important for the key convergence (11.81) . since possible sign changes (if any) 
disappear in the limit as t — > T~ . It seems that nothing prevents existence of some ODEs 
from the family (I2.2p - (I2.6I) . with different nonlinearities, for which the shock profiles can 
change sign for z < 0. 

Proof. As above, we consider the first ODE (12.21) only. We use a shooting argument using 
the 2D bundle of asymptotics (I2.2ip . By scaling (I2.22p . we put C = —1, so, actually, we 
deal with the one-parameter shooting problem with the ID family of orbits satisfying 

(2.24) g(z; D) = -z + D z 3 - ^ z 5 - JjL z 7 + ... , Del. 

It is not hard to check that, besides stabilization to unstable constant equilibria, 

(2.25) g{z) -»• C_ > as z ->• -oo, 
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the ODE (I2.3P admits an unbounded stable behaviour given by 



(2.26) g(z) ~ g*(z) = — ^ z 5 — > +00 as z ^ -00. 

The overall asymptotic bundle about the exact solution g*{z) is obtained by linearization: 
as z — ■> —00, 

(2.27) g(z) = g*(z) + Y(z) => ^r( 5 ) = -±Y'z + ... or z 5 F( 5 ) = 24F'z + ... . 
This is Euler's type homogeneous equation with the characteristic equation 

(2.28) Y(z) = z m m x = (Y 2 (z) = 1) or (m- l)(m- 2)(m- 3)(m - 4) = 24. 

This yields another m 2 = (hence there exists Y 2 (z) = In |z|), m 3 = 5 (not suitable), and 
a proper single complex root with Rem = | < 5 yielding oscillatory Y 3 ^(z). Thus: 

(2.29) as z — > —00, there exists a 4D asymptotic bundle about <7*(z) = — 2;5 - 

Therefore, at z — —00, we are given a 2D bundle of proper solutions (12.171) . as well as 
a 4D fast growing profiles from (12.291) . This determines the strategy of the ID shooting 
via the D-family (12231) : 

(i) obviously, for all D ^ —1, we have that g(z;D) > is monotone decreasing and 
approaches the stable behaviour (I2.26p . (12.291) . and 

(ii) on the contrary, for all D ^> 1, g{z\ D) gets non-monotone and has a zero at some 
finite z = Zq(D) < 0, satisfying z (D) — >• _ as D — >• +00, and eventually approaches the 
bundle in (12.291) . but in an essentially non-monotone way. 

It follows from different and opposite "topologies" of the behaviour announced in (i) 
and (ii) that there exists a constant D such that g(z; D ) does not belong to those two 
sets of orbits (both are open) and hence does not approach g*(z) as z — > —00 at all. This 
is precisely the necessary shock similarity profile. □ 

This ID shooting approach is explained in Figure H] obtained numerically, where 

(2.30) D = 0.069192424.... 

It seems that as D — » Dq, the zero of g(z; D) must disappear at infinity, i.e., 

(2.31) zq(D) -> -00 as D -> D+ , 

and this actually happens as Figure H] shows. Then this would justify the positivity 
(12.231) . Unfortunately, in general (i.e., for similar ODEs with different sufficiently arbitrary 
nonlinearities), this is not true, i.e., cannot be guaranteed by a topological argument. So 
that the actual operator structure of the ODEs should be involved in the study, so, 
theoretically, the positivity is difficult to guarantee in general. Note again that, if the 
shock similarity profile g(z) had a few zeros for z < 0, this would not affect the crucial 
convergence property such as (II. 8p . 
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Shooting shock similarity profile for C=-1 : D =0. 0691 92424... 
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FIGURE 4. Shooting the shock similarity profile g(z) via the family (f!T2"4"j) ; D = 0.069192424... . 

2.5. Self-similar formation of other shocks. NDE-{l,A). Let us first briefly consider 
the last ODE (12. 6p for the fully divergent NDE (11.51) . Similarly, by the same arguments, we 
show that, according to (12.11) . there exist other non-symmetric shocks as non-symmetric 
step-like functions, so that, as t — > CT, 

( C_ > for x < 0, 

(2.32) u-(x,t) -> < C for x = 0, 

{ C+ < for x > 0, 

where C_ ^ —C + and Co 7^ 0. Figure [5] shows a few of such similarity profiles g{z), where 
three of these are strictly positive. The most interesting is the boldface one with 

C_ = 1.4 and C+ = 0, 

which has the finite right-hand interface at z = zq ~ 5, with the expansion 

(2.33) g(z) = (z - z)%{\ + o(l)) ^0" asz^ z . 

It follows that this g(z) < near the interface so the function changes sign there, which 
is also seen in Figure [5] by carefully checking the shape of profiles above the boldface one 
with the finite interface bearing in mind a natural continuous dependence on parameters. 

NDE-(5,0). Consider next the first ODE Q for the fully non-divergent NDE-(5, 0) 
(II. ip . We can again describe formation of shocks (12.321) ; see Figure El The boldface 
profile with C_ = 1.4 and C + = has finite right-hand interface at z = Zq ~ 5, with a 
different expansion 

(2.34) g{z) = ^f{z - zf\ ln(z - z)\{\ + o(l)) -> 0+ as z -> z^. 
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Shock similarity profiles for NDE-(1 ,4)=(0,5): u=- (u u x ). 




2.6. Shock formation for a uniformly dispersive NDE: an example. Here, as a 
key example to be continued, we show shocks for uniform (non-degenerate) NDEs, such 
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Shock similarity profiles for u ( =- ((1+u )u x ) xxxx 
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FIGURE 7. Various shock similarity profiles g(z) satisfying the ODE (j2.36[) . 



as the fully divergent one, 
(2.35) 



u t = -((1 +u 2 )u x ) 3 



where the dispersion coefficient —(1 + u 2 ) of the principal operator is an even function. 
Recall that, for all the previous ones (ll.ip - (ll.5l) . the dispersion coefficient is equal to 
— u and is an odd function of u. Equation (12.351) is non-degenerate and represents a 
"uniformly dispersive" NDE. The ODE for self-similar solutions (12.11) then takes the form 

(2.36) ((l + ^ y) (4) = _i^_ 

The mathematics of such equations is similar to that in Section 127^1 In Figure[?l we present 
a few shock similarity profiles for (I2.36p . Note that both shocks S±(x) are admissible, 
since for the ODE (12. 36ft (and for the NDE (12.351) ). we have, instead of symmetry (12. 8ft . 



(2.37) 



— g(z) is also a solution. 



(3.1 



3. (II) Riemann Problem S + : SIMILARITY RAREFACTION WAVES 
Using the reflection symmetry of all the NDEs (ll.ip - (ll.5l) . 

u i— > — u, 



t ^ -t, 



we conclude that these admit global similarity solutions defined for all t > 0, 



(3.2) 



u + (x, t) = g(z), with z = x/t->. 
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Then g(z) solves the ODEs (I2.2l) - (l2.6p with the opposite terms 

(3.3) ... = \g'z 

on the right-hand side. The conditions (12. 7p also take the opposite form 

(3.4) /(±oo) = ±1. 

Thus, these profiles are obtained from the blow-up ones in ( 12. II) by reflection, i.e., 

(3.5) if g(z) is a shock profile in (12. ip . then g(—z) is a rarefaction one in (13.21) . 

These are sufficiently regular similarity solutions of NDEs that have the necessary initial 
data: by Proposition I2.2( ii). in L\ oc , 

(3.6) u+(x, t) — > S + (x) as t — > + . 

Other profiles g{— z) from shock wave similarity patterns generate further rarefaction 
solutions including those with finite left-hand interfaces. 

4. (Ill) Riemann Problem £>_: TOWARDS 5-ENTROPY TEST 

4.1. Uniform NDEs. In this section, for defmiteness, we consider the fully non-divergent 
NDE (O), 

(4.1) u t = A(u) = -uu xxxxx in Ex (0, T), u(x, 0) = u (x) G C °°(M). 

In order to concentrate on shocks and to avoid difficulties with finite interfaces or transver- 
sal zeros at which u = (these are weak discontinuities via non- uniformity of the PDE), 
we deal with strictly positive solutions satisfying 

(4.2) i < u < C, where C > 1 is a constant. 

Remark: uniformly non-degenerate NDEs. Alternatively, in order to avoid the 
assumptions like (14. 2p . we can consider the uniform equations such as (cf. (12.351) ) 

(4.3) u t = -(1 + u 2 )u XXXXXl 

for which no finite interfaces are available. Of course, (14. 3p admits analogous blow-up 
similarity formation of shocks by (12.11) . In Figure we show a few profiles satisfying 

(4.4) (l+g 2 )g (5) = -\g'z, ^l. 

Recall that, for (14.41) . (12.371) holds, so both S±(x) are admissible and entropy (see below). 

4.2. On uniqueness, continuous dependence, and a priori bounds for smooth 
solutions. Actually, in our 5-entropy construction, we will need just a local semigroup 
of smooth solutions that is continuous is L\ oc . The fact that such results are true for 
fifth-order (or other odd-order NDEs) is easy to illustrated as follows. One can see that, 
since (14. ip is a dispersive equation, which contains no dissipative terms, the uniqueness 
follows as for parabolic equations such as 

Ut = —uu xxxx or u t = uu xxxxxx (in the class < u < Cj). 
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Shook similarity profiles for u=- (1+u )u 

1 r t v ' xxxxx 
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FIGURE 8. Various shock similarity profiles g{z) satisfying the ODE (14. 4 



Thus, we assume that u(x,t) solves (14. Xp with initial data Uq(x) G H 10 (R), satisfies 
(O]) . and is sufficiently smooth, u G L°°([0, T], H 10 (R)), u t G L°°([0, T], H 5 (R)), etc. 
Assuming that v(x, t) is the second smooth solution, we subtract equations to obtain for 
the difference w = u — v the PDE 



(4.5) 



w t = —uw 2 



,w. 



We next divide by u > ^ > and multiply by w in L 2 , so integrating by parts that 
vanishes the dispersive term w xxxxx yields 



(4.6) 



[■ wuh — 1 d_ f i»! i i f u t nn 2 
J u — 2 dt J « f 2 J 7 1 " 



V xxxxx <y^2 



Therefore, using (14. 2 p and the assumed regularity yields 



(4.7) 



1 d_ 

2 dt 



Jm J \ 2 u A u I — 1 J u 



where the derivatives u t (-,t) and v xxxxx (-, t) are from L°°([0, T]). By Gronwall's inequality, 
(14. T[) yields w(t) = 0. Obviously, these estimates can be translated to the continuous 
dependence result in L 2 and hence in L\ oc . 

Other a priori bounds on solutions can be also derived along the lines of computations 
in [H §§ 2,3] that lead to rather technical manipulations. The principal fact is the same 
as seen from (14 .7p : differentiating a times in x equation (14. X p and setting v = D"u yields 
the equations with the same principal part as in (14.51) : 



(4.8) 



'''t xxxxx T 
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Multiplying this by with ( being a cut-off function, and using various interpolation 
inequalities makes it possible to derive necessary a priori bounds and hence to observe 
the corresponding smoothing phenomenon for exponentially decaying initial data. 

4.3. On local semigroup of smooth solutions of uniform NDEs and linear oper- 
ator theory. We recall that local C°°-smoothing phenomena are known for third-order 
linear and fully nonlinear dispersive PDEs; see [7J [SJ [27J E3] and earlier references therein. 
We claim that, having obtained a priori bounds, a smooth local solution can be con- 
structed by the iteration techniques as in [HI § 3] by using a standard scheme of iteration 
of the equivalent integral equation for spatial derivatives. We present further comments 
concerning other approaches to local existence, where we return to integral equations. 

We then need a detailed spectral theory of fifth-order operators such as 

(4.9) P 5 = a{x)D 5 x + b(x)D* + xe(-L,L) (a{x)>±>0), 

with bounded coefficients. This theory can be found in, e.g., Naimark's book [351 Ch. 2]. 
For regular boundary conditions (e.g., for periodic ones that are regular for any order, 
which suits us well), operators (14. 9B admit a discrete spectrum {A^}, where the eigenvalues 
Afc are all simple for all large k. 

It is crucial for further use of eigenfunction expansion techniques that the complete in 
L 2 subset of eigenfunctions {ipk} creates a Riesz basis, i.e., for any / G L 2 , 

(4.10) Eia^)| 2 <oo, where (fM = f f^, 

and, for any {c^} G I 2 (i.e., £ \ck\ 2 < oo), there exists a function / G L 2 such that 

(4.11) (f,ipk) = c k . 

Then there exists a unique set of "adjoint" generalized eigenfunctions {4>l} (attributed 
to the"adjoint" operator P|) being also a Riesz basis that is bi-orthonormal to {ipk}'- 

(4.12) (ipk,^) = 8ki (Kronecker's delta). 
Hence, for any / G L 2 , in the sense of the mean convergence, 

(4.13) / = £c fc ^, with Ck = (f,r k }- 

See further details in [351 § 5]. 

The eigenvalues of ( 14. 9ft have the asymptotics 

(4.14) A fc ~ (±2vrA;i) 5 for all k » 1. 

In particular, it is known that P5 has compact resolvent, which makes it possible to use 
it in the integral representation of the NDEs; cf. [8, § 3], where integral equations are 
used to construct a unique smooth solution of third-order NDEs. 

On the other hand, this means that P5 — al for any a ^> 1 is not a sectorial operator, 
which makes suspicious using advanced theory of analytic semigroups [H [151 El], as is 
natural for even-order parabolic flows; see further discussion below. Analytic smoothing 
effects for higher-order dispersive equations were studied in [51]. Concerning unique 
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continuation and continuous dependence properties for dispersive equations, see [12] and 
references therein, and also [52J for various estimates. 



4.4. Hermitian spectral theory and analytic semigroups. Let us continue to dis- 
cuss related spectral issues for odd-order operators. For the linear dispersion equation 
with constant coefficients ( 12.131) . the Cauchy problem with integrable data uq(x) admits 
the unique solution 

(4.15) u(x,t) = b(x- -,t) *u (-), 

where b(x,t) is the fundamental solution (12.141) . Analyticity of solutions in t (and x) can 
be associated with the rescaled operator 



(4.16) B 5 = -D\ + \zD z + \l in Lj(R), where p(z) = j 



e a N 5/4 , z<0, 
e~ az5/ \ z>0, 



and a > is a sufficiently small constant. Here, B5 in (14.161) is the operator in ( 12.151) that 
generates the rescaled kernel F of the fundamental solution in (12.141) . 
Next, using in (12.131) the same rescaling as in ( 12.141) . we set 

(4.17) u(x,t) = t~^v(y,r), y = x/t$, r = lni, 
to get the rescaled PDE with the operator ( 14.161) . 

(4.18) v T = B 5 v. 
Next, on Taylor expansion of the kernel in (14.151) yields 



(4.19) v(y, t) = J F(y - ze~^) Uo (z) dz = E ^r F(k) (y>~ T ^ I z k u {z) dz, 

(k) 

where the series converges uniformly on compact subsets, defining an analytic solution, 
and also in the mean in L 2 p . According to the eigenfunctions expansion ( 14. 19ft of the 
semigroup, there is a proper definition of the operator ( 14.16!) with a real spectrum and 
eigenfunctions (see details in [T71 § 9], [201 § 8-2]) 

a(B 5 ) = {-|, fc = 0,1,2,...} and fa(y) = F^{y), k > 0. 

The basis of the "adjoint" operator (cf. (12.121) ). in a space with an indefinite metric, 

B% = -Dl-\yD y in Lj.(R), p*(z) = e~ a ^ 5/4 in R, 

has the same point spectrum and eigenfunctions {V^}, which are generalized Hermite 
polynomials. Cf. a full "parabolic" version of such a Hermitian spectral theory in [T3| [TT] . 
This implies that B 5 — al is sectorial for a > (Ao = is simple), and this justifies the 
fact that (I4.15P is an analytic (in t) flow. Let us mention again that analytic smoothing 
effects are well known for higher-order dispersive equations with operators of principal 
type, [EE]. 
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Actually, this also suggests to treat (14 .11) . (14.21) by a classic approach as in Da Prato- 
Grisvard [9J by linearizing about a sufficiently smooth u = u(t ), t > 0, by setting 
u{t) = uq + v(t) giving the linearized equation 

(4.20) v t = A'(u )v + A{u )+g{v), t > t ; v(t ) = 0, 

where g(v) is a quadratic perturbation. Using the good semigroup e A ( UQ )* ( this makes it 
possible to study local regularity properties of the corresponding integral equation 

(4.21) v{t) = j e A '( u °)^ s ) (A(uo) + g(v(s))) ds. 

to 

Note that this smoothing approach demands a fast exponential decay of solutions v(x,t) 
as x — > oo, since one needs that v(-,t) e L 2 p ; cf. [33], where C°°-smoothing for third-order 
NDEs was also established under the exponential decay. Equation (14.211) can be used to 
guarantee local existence of smooth solutions of a wide class of odd-order NDEs. 

Thus, we state the following conclusion to be used later on: 

any sufficiently smooth solution u(x, t) of ( 14. Ill . ( 14. 2ft at t — t , 
can be uniquely extended to some interval t £ (t ,t + v), v > 0. 

4.5. Smooth deformations and ^-entropy test for solutions with shocks. The 

situation dramatically changes if we want to treat solutions with shocks. Namely, it is 
known that even for the NDE-3 (ll.30p . the similarity formation mechanism of shocks 
immediately shows nonunique extensions of solutions after a typical "gradient" catastro- 
phe [21] . Therefore, we do not have a chance to get, in such an easy (or any) manner, 
a uniqueness/entropy result for more complicated NDEs such as (11. 5p by using the 5- 
deformation (evolutionary smoothing) approach. However, we will continue using these 
ideas, turned out to be fruitful, in order to develop a much weaker "5-entropy test" for 
distinguishing some simple shock and rarefaction waves. 

Thus, given a small 5 > and a sufficiently small bounded continuous (and, possibly, 
compactly supported) solution u(x, t) of the Cauchy problem (14. ip . satisfying (14.21) . we 
construct its smooth 5 -deformation, aiming get smoothing in a small neighbourhood of 
bounded shocks, as follows. Note that we deal here with simple shock configurations 
(mainly, with 1-shock structures), and do not aim to cover more general shock geometry, 
which can be very complicated; especially since we do not know all types of simple single- 
point moving shocks. 

(i) We perform a smooth 5-deformation of initial data u (x) by introducing a suitable 
C 1 function Uqs(x) such that 

(4.23) f\u -u os \<5. 

If Mo is already sufficiently smooth, this step must be abandoned (now and always later 
on). By u\s(x, t), we denote the unique local smooth solution of the Cauchy problem with 
data uos, so that, by (14.221) . the continuous function Ui$(x, t) is defined on the maximal 
interval t e [to,ti(5)), where we denote to = and t\(5) = Aig. At this step, we are 
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able to eliminate non-evolution (evolutionary unstable) initially posed shocks, which then 
create corresponding smooth rarefaction waves. 

(ii) At t — A\s, a shock-type discontinuity (or possibly infinitely many shocks) is sup- 
posed to occur, since otherwise we extend the continuous solution by (14. 22ft . so we perform 
another suitable 5-deformation of the "data" uu(x, A 1( j) to get a unique continuous solu- 
tion U2s(x,t) on the maximal interval t G [ti(S), t2(S)), with t2(S) = Ais + A 2< 5, etc. Here 
and in what follows, we always mean a "^-smoothing" performed in a small neighbourhood 
of occurring singularities only as discontinuous shocks. 

We continue in this manner with suitable choices of each (^-deformations of "data" at 
the moments t = tj(5), when ujs{x,t) has a shock, there exists a tk{5) > 1 for some finite 
k = k(5), where k(S) — > +oo as 5 — > 0. It is easy to see that, for bounded solutions, k(S) 
is always finite. A contradiction is obtained by assuming that tj(5) — > i < 1 as j ; — > oo 
for arbitrarily small 5 > meaning a kind of "complete blow-up" that was excluded by 
assumption of smallness of the data. 

This gives a global 5 -deformation in 1 x [0,1] of the solution u(x,t), which is the 
discontinuous orbit denoted by 

(4.24) u 5 (x, t) = {u j5 {x, t) for t e tj(5)), j = 1, 2, .., k(5)}. 

One can see that this ^-deformation construction aims at checking a kind of evolution 
stability of possible shock wave singularities and therefore, to exclude those that are not 
entropy and evolutionarily generate smooth rarefaction waves. 

Finally, by an arbitrary smooth S- deformation, we will mean the function (14.241) con- 
structed by any sufficiently refined finite partition {tj(5)} of [0,1], without reaching a 
shock of SL-type at some or all intermediate points t = tj(S). 

We next say that, given a solution u(x,t), it is stable relative smooth deformations, or 
simply 5-stable (5 eformation- stable) , if for any e > 0, there exists 5 = 5(e) > such that, 
for any finite <5-deformation of u given by (14.241) . 

(4.25) JJ\u-u 5 \<£. 

Recall that (I4.24p is an 5-orbit, and, in general, is not and cannot be aimed to represent 
a fixed solution in the limit 5 —>■ 0; see below. 

4.6. On 5-entropy solutions. Having checked that the local smooth solvability problem 
above is well-posed, we now present the corresponding definition that will be applied to 
particular weak solutions. Recall that the metric of convergence, L\ oc under present 
consideration, for (11.301) was justified by a similarity analysis presented in Proposition 
12.21 For other types of shocks and/or NDEs, the metric may be different. 

Thus, under the given hypotheses, a function u(x,t) is called a 5-entropy solution 
of the Cauchy problem (14.11) . if there exists a sequence of its smooth (^-deformations 
{u &k , k = 1, 2, ...}, where 5k — > 0, which converges in L\ oc to u as k — > oo. 

This is slightly weaker (but equivalent) to the condition of 5-stability. 
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Remark: 5-entropy solution is unique for ID conservation law. Consider, as a 
typical example, fll.7j) for general measurable L 1 -data. The classical Oleinik-Kruzhkov's 
entropy theory for ( 11.7ft defines the unique semigroup of contractions in L 1 (see [50]), i.e., 
for an arbitrary pair of entropy solutions u(-,t) and v(-,t), in the sense of distributions, 



(4.26) ± \\u(t) -v(t)\\ L i < for a.a. t > 0. 

Consider now the above 5-deformation construction of an orbit {u$} in the case, when 
the entropy solution u(x,t) is continuous a.e. for all t > 0, i.e., shocks have zero measure. 
It means that u$(x,t) for t > is smooth and essentially differs from u(x,t) on a set of 
arbitrarily small measure ~ 5 — > 0. Therefore, under these (possibly, non-constructive) 
assumptions, ( 14. 26ft implies that any smooth (^-deformations in L 1 inevitably lead to the 
unique entropy solution of (11.71) as 5 — > 0. In other words, 



(4.27) for Euler's equation (jl.7p . classic entropy solutions = 5-entropy ones. 



Of course, this is just the trivial consequence of the L 1 -contractivity ( 14. 26ft . which, in its 
turn, is induced by the Maximum Principle. It is also worth mentioning that, somehow, 
(14.261) reflects the fact that the conservation laws such as (11.71) admit the direct algebraic 
solution via characteristics. Indeed, the characteristic method guarantees the unique solv- 
ability in the regularity domain, while the "shocks cut off" can be performed at necessary 
points by the corresponding Rankine-Hugoniot relations. Thus, the entropy conditions 
just describe the correct evolution from initially posed singularities (evolutionary, such 
"rarefaction waves" cannot appear by characteristics). 

Therefore, the absence of the Maximum Principle and absence of any characteristic- 
based approaches for higher-order NDEs recall that a result such as (14.271) cannot be 
expected in principle here. The situation is even more terrible: we will show that any 
uniqueness/entropy results for such NDEs fail always and anyway. 

4.7. 5-entropy test and nonexistent uniqueness. Since, for obvious reasons, the 5- 
deformation construction gets rid of non-evolutionary shocks (leading to non-singular 
rarefaction waves), a first consequence of the construction is that it defines the 5-entropy 
test for solutions, which allows one, at least, to distinguish the true simple isolated shocks 
from smooth rarefaction waves. 

In Section we show that it is completely unrealistic to expect from this construction 
something essentially stronger in the direction of uniqueness and/or entropy-like selection 
of proper solutions. Though these expectations correspond well to previous classical PDE 
entropy-like theories, these are excessive for higher-order models, where such a universal 
property is not achievable at all any more. Even proving convergence for a fixed special 5- 
deformation is not easy at all. Thus, for particular cases, we will use the above notions with 
convergence along a subsequence of 5's to classify and distinguish shocks and rarefaction 
waves of simple geometric configurations: 

4.8. First easy conclusions of 5-entropy test. As a first application, we have: 



Proposition 4.1. Shocks of the type S-(x) are 5-entropy for (4-1)- 
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The result follows from the properties of similarity solutions (I2.ip . with —t i— > T t , which, 
by varying the blow-up time T i— > T + 5, can be used as their local smooth 5-deformations 
at any point t > 0. 



Proposition 4.2. Shocks of the type S + (x) are not 5 -entropy for (4-1)- 



Indeed, taking initial data S+(x) and constructing its smooth <5-deformation via the 
self-similar solution (13.21) with shifting t i— > t + 5, we obtain the global 5-deformation 
{u 5 = u + (x, t + 6)}, which goes away from S + . 

Thus, the idea of smooth 5-deformations allows us to distinguish basic 5-entropy and 
non-entropy shocks without any use of mathematical manipulations associated with stan- 
dard entropy inequalities, which, indeed, are illusive for higher-order NDEs; cf. [21J. We 
believe that successful applications of the ^-entropy test can be extended to any configu- 
ration with a finite number of isolated shocks. However, it is completely illusive to think 
that such a simple procedure could be applied to general solutions, especially since the 
uniqueness after singularity formation cannot be achieved in principle, as we show next. 

In other words, the 5-entropy test allows us to prohibit formation of non 5-deformation 
stable shocks of type S+ and proposes a smooth rarefaction wave instead. However, 
this approach cannot detect a unique shock of the opposite geometry since such a 
formation is principally nonunique. 

5. (IV) Nonuniqueness AFTER SHOCK FORMATION 

Here we mainly follow the ideas from [21] applied there to the NDE-3 (ll.30p . so we 
will omit some technical data and present more convincing analytic and numerical results 
concerning the nonuniqueness. For the hard 5D dynamical systems under consideration, 
numerics becomes more and more essential and unavoidable for understanding the nature 
of such nonunique extensions of solutions. Without loss of generality, we always deal with 
the NDE-5 (jl.5p of the fully divergent form. 

5.1. Main strategy towards nonunique continuation: pessimistic conclusions. 

We begin with the study of new shock patterns, which are induced by other (cf. (12.11) ) 
similarity solutions of (11.51) : 

(5.1) u4x,t) = (-trf(y), y = Hb, /3 = ^, where a G (0, i) and 



(5.2) 



<ffT ) -f3f'y + uf = in R_, /(0) = f'(0) = / (4) (0)=0, 
f(y)=C \y\f (l + o(l)) as y^-oo, C > 0. 



In this section, in order to match the key results in [21], in ( 12.1 ft and later on, we change 
the variables {g,z} i— > {f,y}- In the next Section [HI we return to the original notation. 
The anti-symmetry conditions in (I5.2p allow us to extend the solution to the positive 
semi-axis {y > 0} by the reflection —f(—y) to get a global pattern. 

Obviously, the solutions (12.11) . which are suitable for Riemann problems, correspond to 
the simple case a = in (15.11) . It is easy to see that, for positive a, the asymptotics in 
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f)5.2p ensures getting first gradient blow-up at x = as t — ► , as a weak discontinuity, 
where the final time profile remains locally bounded and continuous: 



(5.3) «_0,(T) 



Cq\x\ & for x < 0, 

— 

-Cq\x\ p ioi X > 0, 



where Co > is an arbitrary constant. Note that the standard "gradient catastrophe", 
u x (0, 0~) = — oo, then occurs in the range, which we will deal within, 

(5.4) | < 1 provided that a < \. 

Thus, the wave braking (or "overturning") begins at t — 0, and next we show that it is 
performed again in a self-similar manner and is described by similarity solutions 

(5.5) u + {x,t)=t a F{y), y = %, (3 = i±s, where 



(5.6) 



- (FF') ^ + (3F'y -aF = in R_ , 
F(0) = F >0, F(y) = C \y\f(l + o(l)) as y -oo, 



where the constant Co > is fixed by blow-up data (15.31) . The asymptotic behaviour 
as y — > — oo in (15. 6p guarantees the continuity of the global discontinuous pattern (with 
F(—y) = —F(y)) at the singularity blow-up instant t = 0, so that 

(5.7) m^(x,0~) = m+(x,0 + ) in E. 

Then any suitable couple {/, F} defines a global solution u±(x, t), which is continuous 
at t — 0, and then it is called an extension pair. It was shown in |21j that, for the typical 
NDEs-3, the pair is not uniquely determined and there exist infinitely many shock-type 
extensions of the solution after blow-up at t — 0. We are going to describe a similar 
nonuniqueness phenomenon for the NDEs-5 such as (11.51) . 

It is worth mentioning that, for conservation laws such as (11.71) . such an extension pair 
{/, F} is always unique; see similarity analysis in [21], § 4]. Of course, this is not surprising 
due to existing Oleinik-Kruzhkov's classic uniqueness-entropy theory [371 EI]- Note again 
that any sufficient multiplicity of extension pairs {/, F}, obtained via small micro-scale 
blow-up analysis of the PDEs, would always lead to a principle nonuniqueness, so this 
approach could be referred to as a "uniqueness test" . 

A first immediate consequence of our similarity blow-up/extension analysis is as follows: 

(5.8) in the CP, formation of shocks for the NDE (11.51) can lead to nonuniqueness. 

The second conclusion is more subtle and is based on the fact that, for some initial data 
at t = (i.e., created by single point gradient blow-up as t — > 0~), the whole admitted 
solution set for t > does not contain any "minimal", "maximal", "extremal" in any 
reasonable sense, or any isolated points, which might play a role of a unique "entropy" 
one chosen by introducing a hypothetical entropy inequalities, conditions, or otherwise. If 
this is true for the whole set of such weak solutions of (11.51) with initial data (15.31) . then, 
for the Cauchy problem, 

(5.9) there exists no general "entropy mechanisms" to choose a unique solution. 
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Infinite shock similarity profiles for cc<0, o/p<0 
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FIGURE 9. Infinite shock similarity profiles as solutions of (|5.2p for a < 0: 
a = —0.05 and a = —0.1. 

Actually, overall, (15. 8p and (15.91) show that the problem of uniqueness of weak solutions 
for the NDEs such as (II. 5p cannot be solved in principal. 

On the other hand, in a FBP setting by adding an extra suitable condition on shock 
lines, the problem might be well-posed with a unique solution, though proofs can be very 
difficult. We refer again to a more detailed discussion of these issues for the NDE-3 (11.301) 
in [21]. Though we must admit that, for the NDE-5 (II. 5p . which induces 5D dynamical 
systems for the similarity profiles (and hence 5D phase spaces), those nonuniqueness and 
non-entropy conclusions are more difficult and not that clear as for the NDEs-3, so some 
of their aspects do unavoidably remain questionable and even open. 

Hence, the nonuniqueness in the CP is a non-removable issue of PDE theory for higher- 
order degenerate nonlinear odd-order equations (and possibly not only for those). The 
nonuniqueness of solutions of (II .5p has some pure dimensional natural features, and, more 
precisely, is associated with the dimensions of "good" and "bad" asymptotic bundles of 
orbits in the 5D phase space of the ODE (15.61) . 

5.2. Infinite shock similarity solutions for a < 0. Let us first note that the blow- 
up solutions (15.11) represent an effective way to describe other types of singularities with 
infinite shocks. Namely, assuming that 



we again obtain the same "data" (15.31) but now w_(0,0~) = oo. We do not study in 
any detail such interesting new singularity phenomena and present Figure M showing that 
such infinite shock similarity profiles do exist. For comparison, we indicate the standard 
5_-type profile for a = 0, which coincides with that in Figure [3](d). 

For the NDE-3 such as (11.301) . the infinite shock similarity solutions in the range (15.101) 
were studied in [23l § 4] in sufficient detail. 

26 



(5.10) 



a < and % 



<0, 



5.3. Gradient blow-up similarity solutions. Consider the blow-up ODE problem 
(15. 2p . which is a difficult one, with a 5D phase space. Note that, by invariant scaling 
(12.221) . it can be reduced to a 4th-order ODE with a also even more complicated nonlin- 
ear operator composed from too many polynomial terms, so we do not rely on that and 
work in the original phase space. Therefore, some more delicate issues on, say, uniqueness 
of certain orbits, become very difficult or even remain open, though some more robust 
properties can be detected rigorously. We will also use numerical methods for illustrating 
and even justifying some of our conclusions. As before, for the fifth-order equations such 
as (15.21) . this and further numerical constructions are performed by the MatLab with the 
standard ode45 solver therein. 

Let us describe the necessary properties of orbits {f(y)} we are interested in. Firstly, 
it follows from the conditions in (15.21) that, for y m 0~, 

(5.11) the set of proper orbits is 2D parameterized by fi = f'(0) < and f 3 = f"'(0). 
Secondly and on the other hand, the necessary behaviour at infinity is as follows: 

(5.12) f(y) = C \y\&(l + o(l)) as y - -oo = f), 

where Co > is an arbitrary constant by scaling (12.221) . It is key to derive the whole 4D 
bundle of solutions satisfying (15. 12ft . This is done by the linearization as y — > — oo: 

, K10 , f(y) = fo(y) + Y(y), where f (y) = C (-y)f 

(5.13) 

=► -C ((-y)?Y)^ + /3Y'(-y) + aY + \ (f*(y))® + ... = 0. 

By WKBJ-type asymptotic techniques in ODE theory, solutions of (I5.13P have a standard 
exponential form with the characteristic equation: 

(5.14) Y(y) ~ e^\ 7 = 1 + i (l - |) > 1 =^ C ( 7 a) 4 = (3, 

which has three roots with non-positive real parts, Rea^ < 0, where a\ < is real and 
conjugate a 2j 3 G iR. Hence, we conclude that: 

(5.15) as y — » — oo, the bundle (15.121) is four-dimensional (including Co). 

The behaviour corresponding the bundle (15.151) gives the desired asymptotics. Indeed, 
by ( 15.12ft . we have the gradient blow-up behaviour at a single point: for any fixed x < 0, 
as t — > _ , where y = x/(—t)P — > — oo, uniformly on compact subsets, 

(5.16) u_(z,t) = (-trf(y) = (-t)«C |^p(l + o(l)) ^ C |x|^. 

Let us explain some other crucial properties of the phase space, now meaning "bad 
bundles" of orbits. First, these are the fast growing solutions according to the explicit 
solution 

(5.17) = > o for y«-l. 
Analogously to (15 . 1 3j) . we compute the whole bundle about (15.171) : 

(5.18) f(v) = Mv)+Y(y) =^ JE ^(y 5 Y)^-PY>y + aY + ... = 0. 
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Blow-up f(y): function h (m) for a=1/9 with five real negative roots 



Blow-up f(y): function h (m) for a=1/9, 1/19, and 3/17 (five real roots) 




(a) a = | 




(b) various a 



FIGURE 10. The polynomial h a (m) in (|5.19p for various a G (0, j): five nega- 
tive roots. 



This Euler equation has the following solutions with the characteristic polynomial: 

(5.19) Y(y) = y m =}► h a (m) = (^+i)(^+2)(m+3)(m+4)(m+ 5 ) _ p m + a = Q 

One root m = —5 is obvious that gives the solution (15.171) . It turns our that this algebraic 
equation has precisely five negative real roots for a from the range (15 .4p . as Figure [TOl 
shows. Actually, (b) explains that the graphs are rather slightly dependent on a. Thus: 

(5.20) the bundle about (15.171) is five-dimensional. 

Second, there exists a bundle of positive solutions vanishing at some finite y — *■ y$ < 
with the behaviour (this bundle occurs from both sides, as y — > y^ to be also used) 

(5.21) f 1 (y) = A^\y~^\(l + o(l)) 1 A > 0, 

is 4D, which also can be shown by linearization about (I5.2ip . Indeed, the linearized 
operator contains the leading term 

(5.22) -A 2 (^\y~^\Y)^ + ... = Y(y) ~ \y - y \§, \y-y \%, \y-yo\l, 
which together with the parameter y < yields 

(5.23) the bundle about (15.211) is four- dimensional. 

Thus, (15.111) . (I5.15p . (I5.20p . and (15.231) prescribe key aspects of the 5D phase space we 
are dealing with. To get a global orbit {f(y),y G M_} as a connection of the proper 
bundles (15.1 1 j) and (I5.15p . it is natural to follow the strategy of "shooting from below" by 
avoiding the bundle (I5.2ip . (15.231) . i.e., using the parameters /i )3 in (15.111) . to obtain 

(5.24) y = -oo. 

It is not difficult to see that this profile f(y) will belong to the bundle (15.151) . The proof 
of such a 2D shooting strategy can be done by standard arguments. By scaling (I2.22p . we 
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Typical shooting of blow-up profile f(y) fora=1/9: f Q = f 2 =f 4 =0, ^=-1 




FIGURE 1 1 . The shooting strategy of a blow-up similarity profile f(y) for a = ^, 
with data /(0) = f"(0) = / (4) (0) = and /'(0) = 1; the shooting parameter is 
/'"(0) = 0.0718040128557... . 

always can reduce the problem to a ID shooting (recall that fo = fi = f± = already): 

(5.25) /i = /'(0) = -l and f 3 = f"'(0) is a parameter. 

By the above asymptotic analysis of the 5D phase space, it follows that: 

(I) for / 3 <C — 1 the orbit belongs to the bundle about (15.171) . and 

(II) for / 3 3> 1, the orbit vanishes at finite y along (I5.2ip . 

Hence, by continuous dependence, we obtain a solution f(y) by the min-max principle 
(plus some usual technical details that can be omitted). Before stating the result, for 
convenience, in Figure [TTJ obtained by the ode45 solver, we explain how we are going to 
justify existence of a proper blow-up shock profile f(y); cf. Figure HI 

Thus, we fix the above speculations as follows: 

Proposition 5.1. (i) In the range ( |5.^D , the problem (\5.2\i admits a shock profile f(y). 

We have the following expectation: (ii) f(y) is unique up to scaling (\2.22\i and is positive 
for y < . This remains an open problem that was confirmed numerically. In [5T] , for the 
NDE-3 (11.301) . the phase space is 3D and a full proof is available. 

In fact, this is a rather typical result for higher-order dynamical systems. E.g., we refer 
to a similar and not less complicated study of a 4th-order ODE [25] , where existence and 
uniqueness of a positive solution of the radial bi-harmonic equation with source: 

(5.26) A 2 r u = u p for r = \x\ > 0, w(0) = 1, w'(0) = u'"(0) = 0, u(oo) = 0, 

was proved in the supercritical Sobolev range p > psob = ^zz> N > A. Here, analogously, 
there exists a single shooting parameter being the second derivative at the origin U2 = 
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FIGURE 12. The odd blow-up similarity profiles f(y) in M_ with a = g, yg, and A. 

w"(0); the value u = «(0) = 1 is fixed by a scaling symmetry. Proving uniqueness of such 
a solution in [25J is not easy and leads to essential technicalities, which the attentive reader 
can consult in case of necessity. Fortunately, we are not interested in any uniqueness of 
such kind. Instead of the global behaviour such as (15.171) . the equation (15.261) admits the 
blow-up one governed by the principal operator u^' + ... = u p (u — > +oo). The solutions 
vanishing at finite point otherwise can be treated as in the family (I). 

More numerics by bvp4c. We next use more advanced and enhanced numerical meth- 
ods towards existence (and uniqueness-positivity, see (ii)) of f(y). Figure [121 shows blow- 
up profiles, with /(0) = 0, constructed by a different method (via the solver bvp4c) for 
convenient values a = |, ^, and ^. Note the clear oscillatory behaviour of such patterns 
that is induces by complex roots of the characteristic equation (15.141) . 
Collapse of shocks: "backward nonuniqueness" . This new phenomenon is pre- 
sented in Figure [131 which shows the shooting from y = 0~ for 

(5.27) a = \ =► | = i. 

This again illustrates the actual strategy in proving Proposition (5HJ However, though the 
phase space looks similar, note that here, as an illustration of another important evolution 
phenomenon, we solve the problem with /(0) 7^ 0, so that there exists a non-zero jump 
of U-(x, t) at x = denoted by [■]: 

(5.28) /(0) = / = 10 =^ [u_(0,t)} = 2f (-t) a ^0 as t - 0~. 

Therefore, this similarity solution describes collapse of a shock wave as t — > _ . 

More numerical results of such types are presented in Figures [H and [151 where we 
use other boundary conditions at y — 0. Note that, being extended for y > in the 
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f^(0) = and the shooting parameter is /'(0) = /"'(0). 

anti-symmetric way, by —f(—y), this will give a proper shock wave solution with the nil 
speed of propagation (see the R H condition (15.421) below). 

In a whole, since all these blow-up profiles satisfy the necessary behaviour as y — > 
— oo as indicated in (15.21) . these create as t — > 0~ the same initial data (15.31) . This 
confirms the following phenomenon of "backward nonuniqueness" : initial data (15.51) 
with gradient blow-up at x = can be created by an infinite number (in fact, by a 2D 
subset parameterized, say, by {/o,/i}) of various self-similar solutions (I5.jp . 

Indeed, such a nonuniqueness is directly associated with the fact that, due to (I5.15p . 
the proper asymptotic bundle as y — ► — oo is 3D (for a fixed Co > 0, we have to subtract 
the dimension via the scaling invariance (12.221) ). Therefore, roughly speaking, shooting 
from y — 0~ with 5 parameters /o = f(0), ... , f± — /^(O) allows a 2D (2 = 5 — 3) subset 
of solutions f(y) with shocks at y — 0. A full justification of such a conclusion requires 
a more careful analysis of the phase space including geometry of two "bad" bundles, 
which we do not perform here concentrating on other more important solutions and true 
nonuniqueness phenomena. 

Stationary solutions with a "weak shock" . The ODE in (15.21) and hence the 
PDE (11. 5p admit a number of simple continuous "stationary" solutions. E.g., consider 

(5.29) a = i»f = 2 : f(v) = VWl si g n ?/ and u(x, t) = ±a/R signx. 
Note that these are not weak solutions of the stationary equation 

(5.30) l(u 2 ) xxxxx = in V. 

The classic stationary solution of (15.301) u(x,t) = ±x 2 is smoother at x = 0. 
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FIGURE 14. Blow-up profiles f(y) for a = 1, with /'(0) = -1, /"(0) = 0; 
/(0) 6 [0, 1.5] being a parameter. 



Shock blow-up profiles f(y) fora=1/9, C Q =1, u(x,0")=|x| 1/2 Shock blow-up profiles f(y) fora=1/9, C Q =1, u(x,0")=|xf 




(a) /'(0) = /"(0) = (b) /"(0) = i, /W(0) = 

FIGURE 15. Blow-up profiles f(y) for a = i, for /'(0) = /"(0) = (a) and 
/"(0) = I, /( 4 )(0) = (b); /(0) is a parameter. 
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We will show that such "weak stationary shocks" as in f!5.29j) also lead to nonuniqueness. 

Remark: an exact solution for a critical a. One can see that the quadratic operator 
B(/) = (//') in (15.21) admits the following polynomial invariant subspace: 

W 6 = Span{l, y, y\ y 3 , y\ y 5 } B(W 6 ) C W 6 . 

Restricting the ODE (15. 2p to We yields an algebraic system, which admits an exact 
solution for the following value of the critical a c : 

(5.31) a = a c = g = 0.202381... =*► 3 f(y) = Cy - § y 5 , CgI. 

Since a c > 0, it does not deliver a "saw" -type blow-up profile (having infinite number of 
positive humps) as it used to be for the NDE-3 (11.301) for a c = — ^; see [23, § 4]. 

5.4. Nonuniqueness of similarity extensions beyond blow-up. As in [2T] for the 
NDEs-3, a discontinuous shock wave extension of blow-up similarity solutions (15. lft . (15.21) 
is assumed to be done by using the global ones (15. 5ft . (15.61) . Actually, this leads to watching 
a whole 5D family of solutions parameterized by their Cauchy values at the origin: 

(5.32) F(0) = F > 0, F'{0) =F 1 <0, F"{0) = F 2 , F"'(0) = F 3 , F (4) (0) = F 4 . 

Thus, unlike (15.111) . the proper bundle in (15.321) is 5D. Note that at y — — oo, the solution 
must have the form 

(5.33) F(y) = C \y\&(l + o(l)) as y - -oo (C > 0). 

As above, the 5D phase space for the ODE in (I5.6P has two stable "bad" bundles: 

(I) Positive solutions with "singular extinction" in finite y, where F(y) — > as y — ► < 
0. This is an unavoidable singularity following from the degeneracy of the equations with 
the principal term FF^ leading to the singular potential ~ -p. As in (15.221) . this bundle 
is 4D, and 

(II) Negative solutions with the fast growth (cf. (I5.17P ): 

(5.34) F m (y) = (1 + o(l)) -oo as y -»• -oo. 

The characteristic polynomial is the same as in (I5.19p . so that the bundle is 5D; cf. (I5.20p . 

Both sets of such solutions are open by the standard continuous dependence of solutions 
of ODEs on parameters. The whole bundle of solutions satisfying (15.121) is obtained by 
linearization as y — > — oo in (15. 6ft : 

, KOC , /( 2 /) = F ( 2 /) + r(y), where F (y) = C (-y)? 

(5.35) 

=► -C ((-y)»Y)^ - (3Y'(-y) -aY + \ (F 2 (y))^ + ... = 0. 
The WKBJ method now leads to a different characteristic equation: 

(5.36) Y(y) ~ e<~y)\ T = l + J(l-|)>l =}► C ( 7 a) 4 = -/?, 

so that there exist just two complex conjugate roots with Re < 0, and hence, unlike ( 15. 15ft . 

(5.37) the bundle (15.121) of global orbits {F(y)} is three-dimensional. 
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Nonexistence of global profile F(y) for a=1/9: F Q =1, F ]= 1,F 2 =F 3 =0 Nonexistence of global profile F(y) for a=1/9: F =1, F 1= F 2 =F 3 =F 4 
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FIGURE 17. Unsuccessful ID shooting of F(y) satisfying (|5.6p from y = CP, 
with conditions F(0) = 1, F'(Q) = -1, F"(Q) = F^(0) = 0, and F"'(0) = 
-0.2000223777... being a parameter. 



However, the geometry of the whole phase space and the structure of key asymptotic 
bundles change dramatically in comparison with the blow-up cases, so that the standard 
shooting of positive global profiles F(y) by the ode45 solver yields no encouraging results. 
We refer to Figure [HI which illustrates typical negative results of a standard shooting. 
Figure [T7| looks better and presents shooting a kind of "separatrix" , which however does 
not belong to the necessary family as in (15. 6p . Actually, this means that a ID shooting is 
not possible, and, as we will see, there occurs a more complicated 2D one, i.e., using two 
parameters. 
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y 

FIGURE 18. Global profiles F{y) of <^M> for a = ±, C = 1, F'(0) = 0; F(0) G 
[1, 9] being a parameter. 



Therefore, we now use the bvp4c solver, and this gives the following results for the case 
(15.271) . with C = 1, as usual. Namely, we show that there are two parameters, say, 

(5.38) F = F(0) and F l = F'(0), 

such that, for their arbitrary values from some connected subset in M 2 , including all points 
with F > and F% < 0, the problem (15.61) admits a solution. This is confirmed in Figure 
[Hfor the case F'(0) = and in Figured! for the cases F'(0) = +1 (a) and F'(0) = -1 
(b). Obviously, all these profiles are different and exhibit fast and "non-oscillatory" 
convergence as y — > — oo to the "good" bundle as in (15. 6p with Co = 1. 

Finally, carefully analyzing the dimensions of all the "bad" and "good" asymptotic 
bundles indicated in (i) and (ii) above, plus (15.371) . unlike the result for blow-up profiles 
in Proposition 15.11 we arrive at even stronger nonuniqueness: 



Proposition 5.2. In the range ( |5.^P and any fixed Co > 0, the problem (\5.6§ admits a 



2D family of solutions, which can be parameterized by Fq and F\ . 

Recall again that, for any hope of uniqueness, the extension pair {/, F} must be unique 
(or at least their subset should contain some "minimal" and/or isolated points as proper 
candidates for unique entropy solutions) for any fixed constant Co > 0, which defines 
the "initial data" (15. 3p at the blow-up time t = 0~. This actually happens for the Euler 
equation (\1.7\i ; see [21j § 4], where the similarity analysis is indeed easier and is reduced to 
algebraic manipulations, but not that straightforward anyway even for such a "first-order 
NDE" . 
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5.5. "Initial nonuniqueness" . A new "nonuniqueness" phenomenon is achieved for the 
values of parameters 

(5.39) F(0) = F < and F'(0) = F x > 0. 

Figure l2"Ul a). (b) shows such shock profiles leading to the nonuniqueness, obtained by a 
standard ID shooting via the ode45 solver. Here, two similarity profiles F(y) are obtained 
via distinct types of shooting: relative to the parameter F'(0) = F"'(0) in (a), and relative 
F a (0) in (b). 

The proof of existence of such profiles F is based on the same geometric arguments 
as that of Proposition 12.31 (with the evident change of the geometry of the phase space). 
These two different profiles posed into the similarity solutions (I5.5P show a nonunique way 
to get solutions with initial data (Co = 1 by scaling) at t = + : 

a 

(5.40) uq{x) = |xpsignx in K, 

which already have a gradient blow-up singularity at x — 0. This is another potential 
type of nonuniqueness in the Cauchy problem for (11. 5ft . showing the nonunique way of 
formation of shocks from weak discontinuities, including the stationary ones as in (j5.29p . 

However, bearing in mind Proposition 14.21 saying that the shocks of S^-type are not 
5-entropy (i.e., not stable relative small smooth deformations), one can expect that the 
shocks as in (I5.39|) are also unstable. Indeed, smooth extensions of weak pointwise shocks 
(15.401) via rarefaction self-similar waves given by (15.61) are 5-entropy. In Figure [211 we 
show such a global rarefaction profile F(y) for a — „, which describes smooth collapse of 
the "weak equilibrium" (1 5 . 2 9 [) . One can see that such rarefaction profiles satisfy F(y) = 
~f(y)> where / are the corresponding blow-up ones, as shown in Figure [12] for various a. 

Overall, it seems that the 5-entropy test rules out such an "initial nonuniqueness" with 
data of type S+ as in (I5.40p . where a unique smooth rarefaction extension is available. 
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Shooting global profile F(y) fora=1/9: F.=-l, F =F =-0. 115526..., F 2 =F 4 =0 Global profile F(y) fora=1/9: F =-1, F =0, F =-0. 16648, F 3 =F 4 =0 




(a) F'(Q) = F"'(0) = -0.115526... (b) F'(0) = 0, F"(Q) = -0.16648... 

FIGURE 20. Shooting a proper solution F(y) of (|5.6jl for a = h with data 

F(0) = -1, F 4 = 0, and F'(0) = F"'(0) = -0.115526... (shooting parameter), 

F 2 = (a), and F(0) = -1, ^'(0) = 0, F 2 = -0.16648... (shooting parameter), 
F 3 = F 4 = (b). 



Smooth rarefaction profile F(y) fora=1/9, C =-1 




y 

FIGURE 21. Global rarefaction profile F(y) of <^E) for a = J, C = -1; F(0) = 
F"(0) = F( 4 )(0) = 0. 

On the other hand, for other classes of data of SL-shape (according to Proposition I4.ip . 
such a nonuniqueness can take place; see Section 15.41 

5.6. More on nonuniqueness and well-posedness of FBPs. The nonuniqueness 
(15. 8p in the Cauchy problem (11.51) . (15.31) is as follows: any F(y) yields the self-similar 
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continuation (I5.5p . with the behaviour of the jump at x 



(profiles F{y) as in Figure 



(5.41) 



-[u + (x,t)]\ x=0 = -(u + (0 + ,t)-u + (0~,t)) = 2F t a <0 for t > 0. 



In the similarity ODE representation, this nonuniqueness has a pure geometric-dimensional 
origin associated with the dimension and mutual geometry of the good and bad asymp- 
totic bundles of the 5D phase spaces of both blow-up and global equations. Since these 
shocks are stationary, the corresponding Rankine-Hugoniot (R-H) condition on the speed 
A of the shock propagation: 



(5.42) 



A 



[(uU x )x 



[(/ 2 ) 



2d 



\x=0 



2[/] \y=0 







is valid by anti-symmetry. As usual, (I5.42p is obtained by integration of the equation 
(11.301) in a small neighbourhood of the shock. The R-H condition does not assume any 
novelty and is a corollary of integrating the PDE about the line of discontinuity. 

Moreover, the R-H condition (I5.42p also indicates another origin of nonuniqueness: a 
symmetry breaking. Indeed, the solution for t > is not obliged to be an odd function of 
x, so the self-similar solution (15.51) for x < and x > can be defined using ten different 
parameters {F , i 7 ^}, and the only extra condition one needs is the R-H one: 

(5.43) [(FF')'"}(0) = 0, i.e., F ~ F± + AF{ F^ + 3(F 2 ~) 2 = F+F+ + 4F+F 3 + + 3(F+) 2 - 

This algebraic equations with ten unknowns admit many other solutions rather than the 
obvious anti-symmetric one: 

p- — _p+ p- _ p+ p- — _p+ p~ — p+ and F~ — — F + 

Finally, we note that the uniqueness can be restored by posing specially designed con- 
ditions on moving shocks, which, overall guarantee the unique solvability of the algebraic 
equation in (I5.43P and hence the unique continuation of the solution beyond blow-up. 
This construction is analytically similar to that for the NDEs-3 (11.301) in [21] . 



6. Shocks for an NDE obeying the Cauchy-Kovalevskaya theorem 

In this short section, we touch on the problem of formation of shocks for NDEs that 
are higher-order in time. Instead of studying the PDEs such as (cf. [18, 23J) 



(6.1) 



Uttt 



./• ' •■■xxxi etc. , 



we consider the fifth-order in time NDE (II. lip , which exhibits certain simple and, at the 
same time, exceptional properties. Writing it for W = (u,v,w,g, h) T as 



(6.2) 



u t ■ 
v t -- 
w t 
9t - 



Vx, 
W x , 

■ 9x, 
hx j 
- uu 1 



or Wt = AW X , with the matrix A 



10 

10 

1 

1 
hOOOO 



:!cS 



(11.111) becomes a first-order system with the characteristic equation for eigenvalues 

-A 5 + u = 0. 

Hence, for any u ^ 0, there exist complex roots, so that advanced results on hyperbolic 
systems [TJ [10] cannot be applied. 

6.1. Evolution formation of shocks. For (11.111) . the blow-up similarity solution is 

(6.3) U-(x, t) — g(z), z — x/(—t), where 

(gg')W = (z 5 g')W = 120g'z + 24Qg"z 2 + 120g"'z 3 

(6.4) 

+20g^z 4 + g^z 5 in R, /(too) = ±1. 
Integrating ( 16 Ah four times yields 

(6.5) gg' = z 5 g' + Az + Bz 3 , with constants A = (g'(0)) 2 > 0, B = § g'(0)g"'(0), 
so that the necessary similarity profile g(z) solves the first-order ODE 

(6-6) £ = 

By the phase-plane analysis of (16.61) with A > and B = 0, we easily get the following: 

Proposition 6.1. The problem flff.^p admits a solution g(z) satisfying the anti-symmetry 
conditions ( \2. 9h that is positive for z < 0, monotone decreasing, and is real analytic. 

Actually, involving the second parameter B > yields that there exist infinitely many 
shock similarity profiles. The boldface profile g(z) in Figure [221 (by ( 16.31) . it gives S-(x) 
as t — > 0~~) is non-oscillatory about ±1, with the following algebraic rate of convergence 
to the equilibrium as z — » — oo: 



g{*) 



1 + ^ + - for B = 0, 



.l + f + ••• for B > 0. 
Note that the fundamental solutions of the corresponding linear PDE 
(6.7) 

is also not oscillatory as x — > ±oo. This has the form 

b(x,t) =t 3 F(y), y = x/t, so that b(x, 0) = ... = b m (x, 0) = 0, b U tt(x, 0) = 8(x). 

The linear equation (16. 7h exhibits some features of finite propagation via TWs, since 

u(x, t) = f(x - At) => -A 5 / (5) = / (5) , i.e., A = -1, 

since the profile f(y) disappears from the ODE. This is similar to some canonical equations 
of mathematical physics such as 

Ut = u x (dispersion, A = — 1) and u tt = u xx (wave equation, A = ±1). 

The blow-up solution (16.31) gives in the limit t — > _ the shock S- (x) , and (11.81) holds. 
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)(4) 



10 

z 

FIGURE 22. The shock similarity profile satisfying (|6.4p . 

Since (II .111) has the same symmetry (13. ip as (11.301) . similarity solutions (16.31) . with 
— 1 1 — > t and g(z) i— ► g(—z) according to (13. 5p . also give the rarefaction waves for 5+(x) as 
well as other types of collapse of initial non-entropy discontinuities. 

6.2. Analytic 5-deformations by Cauchy— Kovalevskaya theorem. The great ad- 
vantage of the equation (II .lip is that it is in the normal form, so it obeys the Cauchy- 
Kovalevskaya theorem [521 p. 387]. Hence, for any analytic initial data u(x,0), u t (x,0), 
u u (x,0), u U t(x,0), and u tt tt(x,0), there exists a unique local in time analytic solution 
u(x,t). Thus, ( II .lip generates a local semigroup of analytic solutions, and this makes it 
easier to deal with smooth (^-deformations that are chosen to be analytic. This defines 
a special analytic 5-entropy test for shock/rarefaction waves. On the other hand, such 
nonlinear PDEs can admit other (say, weak) solutions that are not analytic. Actually, 
Proposition 16.11 shows that the shock S-(x) is a 5-entropy solution of (11.111) . which is 
obtained by finite-time blow-up as t — > from the analytic similarity solution (16.31) . 

6.3. On formation of single-point shocks and extension nonuniqueness. Similar 
to the analysis in Section [5j for the model (11.1 ip (and (16.11) ). these assume studying 
extension similarity pairs {/, F} induced by the easy derived analogies of the blow-up 
(I5T2|) and global ([ESD, with 

Q 5+a 
H 5 , 

5D dynamical systems. These are very difficult, so that checking three types (standard, 
backward, and initial) of possible nonuniqueness and non-entropy of such flows with strong 
and weak shocks becomes a hard open problem, though some auxiliary analytic steps 
towards nonuniqueness are doable. Overall, in view of complicated multi-dimensional 
phase spaces involved, we do not have any reason for having a unique continuation after 
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singularity. In other words, for such higher-order NDEs, uniqueness can occur accidentally 
only for very special phase spaces, and hence, at least, is not robust (in a natural ODE- 
PDE sense) anyway. 

7. (V) Problem "Oscillatory Smooth Compactons" of fifth-order NDEs 

We begin with an easier explicit example of nonnegative compactons for a third-order 
NDE. 

7.1. Third-order NDEs: 5-entropy compactons. Compactons as compactly sup- 
ported TW solutions of the K(2,2) equation (11.221) were introduced in 1993, [37], as 

(7.1) u c (x,t) = f c (y), y = x + t f c : / = (f 2 )" + f 2 . 

Integrating yields the following explicit compacton profile: 

( 72 ) f ( v )-fi cos2 {i) for \v\< 2 ^ 

KW-\ o for \y\>2n. 

The corresponding compacton ( 17. 1L ( 17.2ft is G-admissible in the sense of iGel'fand 0(1959) 
[26] §§ 2, 8], and is a ^-entropy solution [18], § 4], i.e., can be constructed by smooth (and 
moreover analytic) approximation via strictly positive solutions of the full third-order 
ODE for f c (y) , 

f = (/T + (f 2 )'- 

Since the PDE is not involved unlike Section I4.5[ the ^-entropy notion coincides with the 
G-admissability. 

It is curious that the same compactly supported blow-up patterns occur in the com- 
bustion problem for the related reaction-diffusion parabolic equation 

(7.3) in = {u 2 ) xx + u 2 . 

Then the standing- wave blow-up (as t — > T~) solution of S-regime leads to the same ODE: 

(7.4) u s (x,t) = (T-t)- 1 /(x) =► f = (f 2 )" + f 2 . 

This yields the Zmitrenko-Kurdyumov blow-up localized solution, which has been known 
since 1975; see more historical details in [23] § 4.2]. 

7.2. Examples of C 3 -smooth nonnegative compacton for higher-order NDEs. 

Such an example was given in [TlT p. 4734]. Following [241 p. 189], we construct this 
explicit solution as follows. The operator F^(u) of the quintic NDE 

(7.5) u t = F 5 (u) = (u 2 ) xxxxx + 25(u 2 ) xxx + lU(u 2 ) x 
is shown to preserve the 5D invariant subspace 

(7.6) W5 = Span{l, cosx, sinx, cos 2x, sin 2x}, 



4.M. Gel'fand, 2.09.1913-5.10.2009. 
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i.e., F 5 (V4 / 5 ) C W 5 . Therefore, (17.51) restricted to the invariant subspace W 5 is a 5D 
dynamical system for the expansion coefficients of the solution 

u(x,t) = d(t) +C 2 (t) cos x + C 3 (t) sin x + C 4 (t) cos2x + C 5 (t) sin2x G W s . 

Solving this yields the explicit compacton TW 

105 cos4 (I) for \y\ ^ n > 



(7.7) u c (x,t) = f c (x + t), where f c (y) 



for \y\ > 7r. 



This C\ solution can be attributed to the Cauchy problem for (17.51) since smooth solutions 
are not oscillatory near interfaces; see a discussion around [2U p. 184]. 

The above invariant subspace analysis applies also to the 7th-order PDE 

(7.8) u t = F 7 (u) = D 7 x (u 2 ) + (3Dl{u 2 ) + ^r(u 2 ) xxx + u(u 2 ) x . 
Here F 7 admits W$ if 

= 25, 7 = 144, and u = 0. 
Moreover [2U p. 190], the only operator F 7 in (17.81) preserving the 7D subspace 

(7.9) W7 = C{1, cos x, sin x, cos2x, sin 2a;, cos3a;, sin3a;} 
is in the following NDE-7: 

(7.10) u t = F 7 (u) = D 7 x (u 2 ) + 77D 5 x (u 2 ) + 1876(u 2 ) xxx + 14400(m 2 ) x . 

This makes it possible to reduce (I7.10p on Wj to a complicated dynamical system. 

7.3. Why nonnegative compactons for fifth-order NDEs are not robust: a 
saddle-saddle homoclinic. Recall that, as usual in dynamical system theory, by ro- 
bustness of trajectories we mean that these are stable with respect to small perturbations 
of the parameters entering the NDE or the corresponding ODEs. In other words, the 
dynamical systems (ODEs) admitting such non-negative "heteroclinic" saddle-like orbits 
0^0 are not structurally stable in a natural sense. This reminds the classic Andronov- 
Pontriagin-Peixoto theorem, where one of the four conditions for the structural stability 
of dynamical systems in M. 2 reads as follows [38l p. 301]: 

(7.11) "(ii) there are no trajectories connecting saddle points... ." 

Actually, nonnegative compactons, such as (17.71) . are special homoclinics of the origin, and 
we will show that the nature of their non-robustness is in the fact that these represent a 
stable-unstable manifold of the origin consisting of a single orbit. Therefore, in consistency 
with (17.111) . the origin is indeed a saddle in M. 4 in the plane {/, /', /", /"'}, obtained after 
integration once; see below. 

In order to illustrate the lack of such a robustness in view of a sole heteroclinic involved, 
consider the NDE (17.51) . for which, substituting the TW solution, on integration, we obtain 
the following ODE: 

(7.12) u c (x,t) = f c (x + t) => / c : 2/ = (f) (4) + -, 
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where we omit the lower-order terms as / —>■ 0. Looking for the compacton profile / > 0, 
we set f 2 = F to get 



(7.13) F (4) = 2VF + ... for y > 0, F'(0) = F"'(0) = 0. 

As usual, we look for a symmetric F(y) by putting two symmetry conditions at the origin. 

Let y = yo > be the interface point of F(y). Then, looking for the expansion as 
y — > y^ in the form 

(7.14) F(y) = ^(y -yf + e(y), with e(y) = o((y - y) 8 ), 
we obtain Euler's equation for the perturbation e(y), 

(7.15) _i_( 2/0 _ 2/ )4 £ (4)_ £ = . 

Hence, e(j/) = (y _ 2/)"\ with the characteristic equation 

(7.16) m(m - l)(m - 2)(m - 3) - 840 = m l = -4, m 2;3 = 3±i ^ TTT , m 4 = 7. 

Hence, Rem; < 8, and, in other words, (17.151) does not admit any nontrivial solution 
satisfying the condition in (17.141) ; see further comments in [241, p. 142]. In fact, it is easy 
to see that (17.151) with e = is the unique positive smooth solution of F^ = 2yF. Thus, 

(7.17) the asymptotic bundle of solutions (I7.14p is ID, 

where the only parameter is the position of the interface yo > 0. 

Obviously, as a typical property, this ID bundle is not sufficient to satisfy (by shooting) 
two conditions at the origin in (17. 13ft . so such TW profiles F(y) > are nonexistent for 
almost all NDEs like that. In other words, the condition of positivity of the solution, 

(7.18) to look for a nontrivial solution F > for the ODE in (17. 131) 

creates a free-boundary "obstacle" problem that, in general, is inconsistent. Skipping 
the obstacle condition (17.181) will return such ODEs (or elliptic equations), with a special 
extension, into the consistent variety, as we will illustrate below. 

Thus, nonnegative TW compactons are not generic (robust) solutions of (2m + l)th- 
order quadratic NDEs with m = 2, and also for larger m's, where some kind of (17.171) . as 
a "dimensional defect" (the bundle dimension is smaller than the number of conditions 
at y — to shoot), remains valid. 

7.4. Nonnegative compactons are robust for third-order NDEs only. The third- 
order case m = 1, i.e., NDEs such as (11.221) . is the only one where propagation of pertur- 
bations via nonnegative TW compactons is structurally stable, i.e., with respect to small 
perturbation of the parameters (and nonlinearities) of equations. Mathematically speak- 
ing, then the ID bundle in (17.171) perfectly matches with the single symmetry condition 
at the origin, 

F" = 2VF + ... and F'(0) = 0. 
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Oscillatory compacton TW profiles for u =-(uu ) +uu , X=-1 

3 r r t v x'xxxx X 



F(y) 4 



. F< 4 »=F-2|F| 1/2 sign F , 


— A./ — — — i_ — — — — — — 







-10 -5 5 10 



y 

FIGURE 23. First two compacton TW profiles F(y) satisfying the ODE in ([7719]) . 

7.5. Compactons of changing sign are robust and 5-entropy for the NDEs— 5. 

As a typical example, we consider the perturbed version (11.121) of the NDE-(1,4) (11. 5p . 
As we have mentioned, this is is written for solutions of changing sign, since nonnegative 
compactons do not exist in general. Looking for the TW compacton (I7.12p yields the 
ODE 

(7.19) f = -\ (\f\fp + \ \f\f => =F- 2\F\-lF for F = \f\f. 

Such ODEs with non-Lipschitz nonlinearities are known to admit countable sets of com- 
pactly supported solutions, which are studied by a combination of Lusternik-Schnirel'man 
and Pohozaev's fibering theory; see [22] . 

In Figure [23J, we present the first TW compacton patterns (the boldface line) and the 
second one that is essentially non-monotone. These look like standard compacton profiles 
but careful analysis of the behaviour near the finite interface at y = ya shows that F(y) 
changes sign infinitely many times according to the asymptotics 

(7.20) F(y) = (y -y) 8 lv(s + s ) + o(l)], s = ]n(y - y) as y ^ y . 

Here, the oscillatory component (p(s) is a periodic solution of a certain nonlinear ODE 
and so is an arbitrary phase shift; see [2j, § 4.3] and [221 § 4] for further details. Thus, 
unlike fTTTTTj) . 

(7.21) the asymptotic bundle of solutions (17.201) is 2D (parameters are y and s ), 

and exhibits some features of a "nonlinear focus" (not a saddle as above) on some manifold. 
Hence, this is enough to match also two symmetry boundary conditions given in (I7.13p . 
Such a robust solvability is confirmed by variational techniques that apply to rather 
arbitrary equations such as in (17. 19[) with similar singular non-Lipschitz nonlinearities. 
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Let us also note that such oscillatory compactons are also 5-entropy in the sense that 
can be approximated by analytic TW solutions of the same ODE but having finite number 
of zeros (i.e., admit smooth analytic (^-deformations); see [HI § 8.3] for a related NDE-5. 

Regardless the existence of such sufficiently smooth compacton solutions, it is worth 
recalling again that, for the NDE (I1.12p . as well as f!2.35|) and (14.31) . both containing 
monotone nonlinearities, the generic behaviour, for other initial data, can include forma- 
tion of shocks in finite time, with the local similarity mechanism as in Section [2] and in 
Section 15.11 representing more generic single point shock pattern formations. 



8. Final conclusions 

The fifth-order nonlinear dispersion equations (ll.ll) - (ll.5p (NDEs-5), which are asso- 
ciated with a number of important applications, are considered. The main achieved 
properties of such ID degenerate nonlinear PDEs are as follows: 

(i) Section [21 these NDEs admit blow-up self-similar formation from smooth solutions 
of shock waves of a specific oscillatory structure, which correspond to initial data S_ (x) = 
— signx, i.e., the same as for the first-order conservation law (11.30P ; 

(ii) Section [31 as customary, self-similar rarefaction waves, which get smooth for any 
t > 0, are created by the reversed data S+(x) = sign re; 

(iii) Unlike the classical theory of first-order conservation laws developed in the 1950s- 
60s, the entropy-type techniques are no longer applied for distinguishing general proper 
(unique) solutions. A 5 -deformation test via smoothing the solutions is developed in 
Section HI which is able to separate shocks and rarefaction waves for particular classes 
of initial data ~ S±(x) with a simple geometry of initial shocks; 

(iv) Section [51 by studying more general self-similar solutions of NDEs, it was shown 
that uniqueness of the solutions after formation of a shock is principally impossible. 
Namely, there exist single point gradient blow-up similarity solutions, which admit an 
infinite number of self-similar extensions beyond. This 2D set of shock wave extensions 
after singularity does not have any distinguished solutions (say, maximal, minimal, iso- 
lated, etc.). This also suggests that any entropy-like mechanisms for a unique continuation 
do not exist either. However, using a proper free-boundary setting, i.e., posing special con- 
ditions on shocks, can restore uniqueness; and 

(v) Section [71 nonnegative compacton solutions for some NDEs-5 are shown to be 
non-robust (not "structurally stable"), i.e., these disappear after a. a. arbitrarily small 
perturbations of the parameters (nonlinearities) of the equations. However, oscillatory 
compactons of changing sign near finite interface are shown to exist and to be robust. 
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